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FINAL  REPORT 

Most  of  my  work  concerns  H°°  control  but  focuses  on  substantially  different  parts  of  the  subject,  namely, 
nonlinear  systems,  optimization  theory  and  algorithms  for  frequency  domain  design,  and  computer  algebra  tailored 
to  systems  and  control  research. 

Nonlinear  systems 

The  modern  approach  to  worst  case  design  in  the  frequency  domain  arose  from  studies  of  amplifier  design  the 
“dual”  problem  of  making  a  circuit  dissipative  using  feedback.  For  linear  systems  key  cases  of  this  were  solved  in 
1965  (SISO)  by  Youla  and  Saito  and  (MIMO)  in  1976  by  Helton.  In  the  early  80’s  Zames  and  Francis  formulated 
H°°  control  and  solved  the  math  problem  by  drawing  on  the  earlier  solutions  to  this  circuits  problem.  In  the 
beginning  the  subject  of  H°°  control  evolved  quickly  in  significant  part  because  key  math  problems  were  already 
reasonably  understood  by  operator  theorists.  I  participated  in  this  earlier  work  (e.g.  solved  the  MIMO  H°° 
control  problem  with  Zames  and  Francis,  also  Pearson  and  Chang)  but  at  the  same  time  begin  pushing  in  new 
directions:  nonlinear  plants  and  an  H°°  approach  to  classical  control. 

Of  the  various  solutions  to  CTRL  one  which  is  easy  to  implement  and  numerically  sound  is  the  Doyle-Glover- 
Kargonekar-Francis  DGKF  two  Riccati  equation  solution.  Consequently  extending  this  to  nonlinear  plants  is  of 
considerable  importance.  Just  prior  to  the  contract  period  there  has  been  considerable  progress  by  Isidori  and 
coworkers  and  by  our  group  (Ball  Helton  'Walker  Zhan).  Isidori  et  al  find  local  sufficient  conditions  and  compute 
(with  Krener’s  software)  power  series  solutions  to  model  problems.  All  of  these  approaches  assume  something 
like  the  dimension  of  the  compensator’s  state-space  equals  that  of  the  controller  state-space. 

Evaluating  performance  of  piecewise  linear  systems  is  an  area  where  we  made  progress.  We  took  a  typical 
architecture  (a  la  Campo-Morari)  for  a  system  with  saturation  and  extracted  one  of  the  key  computational 
difficulties.  These  systems  are  piecewise  linear  and  continuous.  Work  with  Ball  showed  that  a  key  object  for  a 
dissipative  system,  called  a  storage  function,  must  be  continuous.  We  then  made  a  natural  compromise.  The 
continuity  of  the  storage  function  forces  constraints  which  make  analyzing  such  systems  not  a  Linear  Matrix 
Inequality.  We  found  a  sequence  of  steps  which  extracted  the  non  LMI  part  and  allowed  one  to  solve  the  problem 
of  determining  performance  of  such  systems  by  doing  first  an  LMI  check,  then  a  side  test  then  an  LMI,  etc. 

General  area  of  nonlinear  H°°  control  In  the  general  area  of  nonlinear  H°°  control  we  settled  some  basic 
theoretical  issues.  James  and  Baras  have  necessary  and  sufficient  conditions  on  the  H°°  control  problem.  James 
and  I  have  extended  the  basics  of  this.  Under  a  saddle  point  assumption  these  reduce  to  cases  also  studied  by  van 
der  Schaft  and  Basar.  Krener  has  results  of  a  similar  tone.  Vityaev  and  I  gave  the  first  theory  on  what  PDE’s 
arise  when  this  saddle  point  structure  fails. 

This  theory  converts  the  problem  of  doing  H°°  control  for  a  nonlinear  system  to  solving  two  particular  PDEs. 
One  PDE  which  computes  optimal  feedback  can  be  solved  off  line.  One  which  gives  the  dynamics  of  the  controller 
must  be  solved  on  line.  Unfortunately,  these  are  PDEs  on  the  state  space  of  the  original  plant  (often  a  high 
dimensional  space)  so  numerical  solution  faces  what  is  called  the  curse  of  dimensionality. 

Beating  the  online  curse  of  dimensionality  My  work  with  James  now  indicates  that,  for  the  mixed  sensitivity 
problem,  the  controller  dynamics  in  practice  might  not  suffer  prohibitively  from  the  curse  of  dimensionality.  In 
this  case  the  biggest  parts  of  the  computation  can  be  done  off  line. 

The  main  observations  leading  to  this  optimism  (for  the  mixed  sensitivity  problem)  are: 

1.  The  controller  PDE  has  some  highly  singular  solutions  pt  which  equal  -oo  off  of  a  small  manifold  Mt.  This 
reduces  the  computational  burden  to  handling  functions  on  Mt. 
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2.  For  linear  plants,  the  DGKF  solutions  can  be  put  in  various  coordinates  while  for  nonlinear  plants  many 
fewer  coordinate  changes  are  possible.  If  one  uses  the  coordinates  which  are  natural  to  the.  nonlinear  theory 
in  a  linear  H°°  situation,  then  one  gets  exactly  the  singular  solutions  in  (1).  Indeed  (1)  gives  the  “central 
controller”  which  solves  the  linear  H°°  problem. 

3.  There  is  a  solid  theory  in  the  nonlinear  case  when  one  uses  smooth  rather  than  singular  solutions  to  the 
controller  PDE. 

4.  What  is  needed  is  to  extend  this  smooth  function  theory  to  singular  functions.  One  thing  we  do  know  is  that 
smooth  solutions  asymptotically  converge  to  a  singular  solution  Poo,  supported  on  the  antistable  manifold 
Mantis  of  the  plant.  This  shows  that  if  a  smooth  solution  to  the  control  problem  exists,  then  Poo  exists  and 
has  some  good  properties.  Now  poo  is  the  most  natural  initialization  for  the  controller  PDE  and  it  produces 
our  singular  controller  in  (1).  Hopefully,  this  controller  solves  the  H°°  problem,  but  this  is  far  from  being 
proved  without  very  strong  assumptions. 

In  conclusion,  once  a  solution  to  the  state  feedback  control  problem  for  a  plant  whose  antistable  dimension  is 
0,  1,  or  maybe  2,  there  are  now  in  principle  formulas  one  could  try  for  solving  the  measurement  feedback  part 
of  the  control  problem.  (The  state  feedback  PDE  remains  oppressive.)  Previously  there  was  a  reasonable  theory 
(as  in  item  3)  but  no  formulas.  Now  there  are  formulas  but  not  much  theory. 

Current  work  with  James  on  theory  gives  (probably  too  conservative)  conditions  along  the  lines  required  in 
(4)  above  on  when  the  singular  controllers  solve  the  control  problem. 

Parameterization  of  all  H°°  controllers  Parameterization  of  all  linear  H°°  controllers  for  a  system  is  equiv¬ 
alent  to  J  inner/outer  factorization  of  the  system.  The  control  and  factorization  problem  for  stable  nonlinear 
systems  was  reduced  to  a  Hamiltion-Jacobi-Bellman-Issacs  (HJBI)  equation  by  Ball  and  Helton  (published  in 
1992).  This  left  the  unstable  case  open. 

Helton  and  James 

1.  gave  formal  equations  along  [JB]  lines  for  J  inner/ outer  factors  and  prove  properties  which  make  their 
formulas  look  very  promising, 

2.  tune  the  [JB]  solution  a  bit  to  correct  for  an  oversight  and 

3.  tighten  the  necessary  [JB]  conditions  to  come  much  closer  to  a  necessary  and  sufficient  theory. 

Optimization  over  H°° 

Much  of  my  effort  goes  to  studying  a  basic  question  of  worst  case  frequency  domain  design  where  stability  of 
the  system  is  the  key  constraint.  This  is  the  H°°  optimization  problem  which  is  crucial  in  several  branches  of 
engineering. 

The  fundamental  H°°  problem  of  control.  First  we  state  the  core  mathematics  problem  graphically.  At 
each  frequency  w  we  are  given  a  set  Su(c)  C  CN,  called  the  specification  set.  The  objective  is  to  find  a  function 
T  with  no  poles  in  the  R.H.P.  so  that  each  T(ju)  belongs  to  «Sw(c).  In  fact  there  is  a  simple  picture  to  think  of 
in  connection  with  a  design 


A“7A 
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T(jw) 


w  =  wt  w  -  wb 


Figure  1 

Typically  there  is  a  nested  family  of  target  sets  <Sw(e)  parameterized  by  a  performance  level  c.  (The  smaller 
the  sets  the  better  the  performance.)  For  the  optimal  c  a  solution  T  exists  but  no  solution  exists  for  tighter  specs. 

The  Horowitz  templates  of  control  can  be  transformed  into  this  type  of  picture.  When  each  <Sw(c)  is  a  “disk” 
this  problem  is  solved  by  transformations  of  “classical  pure”  mathematics  done  in  the  late  1970’s  by  Helton. 
Many  different  solutions  to  this  problem  in  many  different  coordinates  were  worked  out  by  engineers  in  the  last 
15  years  since  it  is  the  subject  of  H°°  control.  Competing  constraints  and  plant  uncertainty  lead  immediately  to 
spec  sets  which  are  not  disks. 

The  graphical  problem  of  Figure  1  can  be  formulated  analytically  in  terms  of  a  performance  function  F  as 
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•  (OPT)  Given  a  positive  valued  function  r  on  R  x  CN  (which  is  a  performance  measure),  find  7*  >  0  and 
/*  in  An  which  solve 

7*=  inf  sup  T(u,f(ju))  . 

f  EAn  oj 

and  this  of  course  is  what  one  puts  in  a  computer.  Collaborators  and  I  have  a  very  broad  based  attack  on  the 
problem  which  addresses  most  aspects  of  it. 

From  qualitative  theory  to  numerical  algorithms  and  diagnostics  While  little  was  known  about  this 
problem  10  years  ago  there  has  been  a  lot  of  progress,  and  now  we  have  a  substantial  amount  of  theory.  We  shall 
not  sketch  all  that  is  known  about  OPT  but  emphasize  that  one  of  the  most  practical  results  on  an  optimization 
problem  is  characterization  of  the  optimum,  since  this  is  the  basis  for  numerics.  We  have,  in  the  last  few  years, 
managed  to  push  from  high  level  theory  to  very  effective  computer  algorithms. 

Time  domain  constraints  Merino,  Walker  and  I  were  able  to  add  time  domain  constraints  to  OPT  and  obtain 
optimality  conditions  extending  those  we  already  had  for  the  OPT  problem.  Our  result  is  easier  to  state  on  the 
unit  disk  A  and  the  unit  circle  T  rather  than  on  the  R.H.P.  and  the  jui- axis.  Also  we  state  it  only  for  the  N=2 
MIMO  case. 

We  consider  a  constrained  optimization  problem,  named  Constr-OPT,  where  the  minimization  is  done  over 
analytic  functions  (/i,/2)  that  satisfy  a  given  set  of  constraints: 

r27T  /*27T 

/  fiGijdO  4-  /  foG2,id6  >  0,  l  —  \  ,  ...,n 

Jo  Jo 

where  the  functions  Gij  are  analytic.  Roughly  the  optimality  condition  for  solutions  to  Constr-OPT  is 

Result  1  Given  T  a  smooth  function  and  the  constraints  above  and  a  smooth  function  T*  in  H°°  satisfying 
a(e10)  =  (et0 ,T*(e10))  is  never  0  on  T.  Necessary  and  sufficient  conditions  for  T*  to  be  a  local  solution  to 
Constr-OPT  are 

IT(et0  ,T*(et0))  is  constant  in  el6 . 

II  There  exist  Fi  and  F2  analytic  on  the  disk ,  A  a  positive  function  on  the  circle ,  and  nonnegative  constants 

K\ , . . . ,  Kn  such  that  for  all  e 10  €  T, 

BT 

(e",  T*{eie))  =  X(eie)  (eieF1(ei9)  +  KlG^  +  . . .  +  nnG^) 

flP 

0^  (eie,  T*{eie))  =  \{eiB)  (ei9F2(ei9)  +  Kl&J  +...  +  n„G ^) 

III  A  postivity  condition  on  second  derivatives  of  T. 

When  there  are  no  time  domain  constraints  this  gives  Kj  =  0  which  is  our  workhorse  result  on  OPT.  Our  analysis 
shows  that  Result  1  meshes  well  the  basic  OPT  result  for  the  purpose  of  constructing  computer  algorithms.  We 
have  worked  out  such  algorithms  and  have  begun  testing. 

Of  independent  interest  is  that  all  of  this  represents  a  new  connection  between  engineering  and  an  existing 
branch  of  the  mathematics  area  Several  Complex  Variables. 

Multiple  performance  objectives  In  H°°  control  one  typically  optimizes  a  supremum  norm  type  of  perfor¬ 
mance  function.  It  has  been  known  for  many  years  that  at  optimum  this  performance  function  is  frequency 
independent  (i.e.  flat). 

We  now  consider  two  competing  performances  Ti  and  T2  which  produces  the  2-OPT  problem. 

Definition  A  function  f *  e  Hff  is  called  a  Pareto  optimum  for  Ti,]^  if  for  each  f  e  Hff,  we  have 

supTi^,/)  >  sup  Pi  (o;,/*)  or  supT2(u;,/)  >  supr2(u;,/*) . 

LJ  UJ  U)  UJ 

In  other  words,  we  cannot  improve  one  of  the  competing  performances  without  degrading  others. 

Andrei  Vityaev  and  I  showed  that  if  the  performance  functions  Ti ,  T2  satisfy  certain  strong  assumptions  and  if 
there  are  N  designable  subsystems  (/x,  ...,/jv)  and  l  performance  measures  with  l  <  N,  then  at  a  nondegenerate 
Pareto  optimum  (/*, ...,  y)  every  performance  is  flat: 

Ti  {ei6,f*(uj))  =  const , . . . ,  Vt(eie,  /»)  =  const . 
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Besides  flatness  there  are  other  “gradient  alignment”  conditions  which  must  hold  at  an  optimum.  Thus  we  have 
the  precise  “first  derivative”  test  for  the  most  natural  class  of  H°°  Pareto  optima. 

MIMO  performance  measures 

Most  recent  work  with  Merino  and  Walker  has  been  in  extending  our  study  of  OPT  to  the  situation  where 
the  performance  T  is  of  the  form 

r(eie,f(eie))  =  \\T(eie,f(eie))\\nxn 

where  T(etG,  f)  is  a  a  smooth  self-adjoint  n  x  n  matrix  valued  function.  This  representation  is  general  enough  to 
cover  practically  all  situations  that  arise  in  applications,  while  at  the  same  time  allowing  the  non-smoothness  of 
OPT  to  be  treated  as  a  feature  of  the  matrix  norm  only.  This  leads  to  very  effective  analysis  of  the  problem  and 
to  algorithms  for  solving  it.  A  key  result  we  obtained  is  the  characterization  of  local  solutions  to  OPT  in  this 
case: 

Result  2  Suppose  that  T  is  matrix  valued .  Under  (mild)  hypotheses ,  if  /*  is  a  local  solution  to  OPT  for  the 
performance  ||r||nxn, 

then  there  exists  a  self  adjoint  n  x  n  matrix  valued  function  $  =  such  that 


(7*1  -r(.,  /*))».  =o 

££ (■,/.)#.)]  =0,  1  <£<N 

f  =1 

>0 

7*7 -r(-,/*)  >0 


Our  characterization  of  local  solutions  leads  immediately  to  algorithms  and  optimality  diagnostics.  The  algo¬ 
rithms,  cases  of  which  we  have  been  working  on  for  about  five  years,  are  of  primal-dual  semidefinite  programming 
(SDP)  type,  so  we  are  merging  our  H°°  optimization  theory  with  SDP  in  i?n. 

We  have  found  that  in  some  cases  it  is  possible  to  have  second-order  convergence  rate,  despite  the  fact  that 
solutions  are  not  unique.  Numerical  analysis  of  the  algorithms  is  in  progress. 

Also,  we  are  exploring  with  G.  Balas  the  efficacy  of  placing  diagnostics  derived  from  our  theory  in  the  Matlab 
package  //-tools.  Such  diagnostics  tell  a  practitioner  how  good  an  approximate  answer  to  a  //-synthesis  problem 
is,  and  aids  further  code  development. 

Computer  algebra  for  systems  research 

If  one  reads  a  typical  article  on  A,B,C,D  systems  in  the  control  transactions  one  finds  that  most  of  the  algebra 
involved  is  noncommutative.  Thus  for  symbolic  computing  to  have  much  impact  on  linear  systems  research  one 
needs  a  program  which  will  perform  noncommuting  operations.  Mathematica,  Macsyma  and  Maple  (the  3  M’s)  do 
not.  We  have  a  package,  NC Algebra,  which  runs  under  Mathematica  and  does  the  basic  operations,  block  matrix 
manipulations  and  other  things.  The  package  might  be  seen  as  a  competitor  to  a  yellow  pad.  Like  Mathematica 
the  emphasis  is  on  interaction  with  the  program  and  flexibility. 

Mins  and  maxes  of  Hamiltonians  Originally  we  wrote  the  package  to  do  linear  H°°  control  research. 
In  particular,  the  main  object  in  studying  CTRL  is  an  energy  balance  (game  theoretic)  Hamiltonian.  One 
must  compute  critical  points  (maxes  or  mins)  of  this  in  W,  a,  6,  c  in  various  orders  which  is  a  routine  but  tedious 
process.  Also  any  variation  on  the  problem  produces  a  new  Hamiltonian  and  requires  another  tedious  computation. 
NCAlgebra  automates  this.  For  example,  if  our  Hamiltonian,  labeled  Ham,  is  quadratic  in  W  and  c,  then 

critW  =  Crit[Ham,W]\  HnoW  —  Ham/ .critW; 
critWc  =  Crit[HnoW ,  cj;  HnoW  c  —  HnoW /  .critW c; 

finds  the  critical  point  of  Ham  in  W  then  in  c  and  and  evaluates  Ham  at  these  critical  points. 

Our  research  focuses  on  what  types  of  “intelligence”  to  put  in  the  package. 

Simplification  of  messy  formulas  We  are  beginning  to  add  serious  automatic  simplification  commands 
to  NCAlgebra.  Stankus,  Wavrik,  and  I  are  now  doing  research  in  computer  simplification  for  A,  B,  C,  D  type 
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linear  systems,  in  a  highly  noncommutative  setting.  The  objective  is  in  each  particular  situation  to  find  a  list  of 
simplifying  rules.  A  complete  list  of  rules  (called  a  Grobner  basis  GB)  has  the  property  that  if  it  is  applied  to  an 
expression  until  nothing  changes  then  the  expression  is  as  simple  as  possible  in  a  certain  sense.  Recently,  Wavrik 
and  I  obtained 


For  the  formulas  which  occur  in  studying  energy  conserving  (lossless)  systems,  the  GB  while  infinite 
can  be  summarized  as  a  list  of  32  rules  some  of  which  depend  on  an  integer  parameter .  It  is  a  powerful 
tool  for  studying  a  particular  class  of  systems .  The  list  was  discovered  last  year  and  actually  proved 
(with  Stankus)  to  be  a  GB  very  recently. 

A  subset  of  these  rules  is  now  in  a  function  NCSimplifyRational[  expression]  inside  our  NCAlgebra  package. 
They  are  very  effective  on  a  limited  class  of  expressions  but  even  that  makes  them  very  useful. 

Discovering  formulas  Current  work  developes  something,  we  call  a  strategy;  it  is  in  a  primitive  stage.  These 
are  methods  for  ”  discovering”  algebraic  theorems  and  formulas  semiautomatically. 

This,  like  the  simplification  methods,  is  based  on  what  is  called  a  noncommutative  Grobner  Basis  Algorithm 
(GBA).  The  GBA  has  the  effect  of  systematically  eliminating  unknowns  so  as  to  put  a  system  of  polynomial 
equations  into  “triangular  form”.  The  commutative  case  of  the  GBA  is  the  core  of  the  “Solve”  commands  of  the 
3  M’s  and  it  is  used  in  many  fields. 

The  input  to  a  GBA  is  (1)  a  list  of  knowns,  (2)  a  list  of  unknowns  (together  with  priorities  for  eliminating 
them)  and  (3)  a  collection  of  equations  in  these  knowns  and  unknowns.  A  strategy  is  :  run  the  GBA,  sort  the 
output  into  equations  involving  only  one  unknown  (say  one  contains  only  aq),  the  user  must  now  make  a  decision 
about  equations  in  aq  (e.g.,  this  is  a  Riccati  so  I  shall  not  try  to  simplify  it,  but  leave  it  for  Matlab).  Now  the 
user  declares  the  unknown  aq  to  be  known  and  runs  the  GBA  again.  Sometimes  one  needs  a  2-strategy  in  that 
the  key  is  equations  in  2  unknowns.  The  point  is  to  isolate  and  to  minimize  what  the  user  must  do.  We  organize 
strategies  via  a  “spreadsheet”  for  discovering  theorems. 

We  are  under  the  impression  that  many  theorems  in  engineering  systems  theory  are  of  this  type.  At  the 
beginning  of  “discovering”  a  theorem,  a  problem  is  presented  as  a  large  system  of  matrix  equations.  Often  when 
viewing  the  output  of  the  GB  algorithm,  one  can  see  what  additional  hypotheses  should  be  added  to  produce  a 
useful  theorem  and  what  the  relevant  matrix  quantities  are.  Our  efforts  are  in  a  primitive  stage  and  the  brevity 
of  this  exposition  suppresses  some  of  the  advantages  and  some  of  the  difficulties. 

Example  Suppose  we  are  given  a  collection  of  equations  involving  matricies.  For  example, 


(FAC) 


Am\  —  mia  —  7712/c  =  0 
Am  2  —  771 2  G  =  0 
B  —  mi  6  —  7712/  =  0 
—c  +  Cmi  =0 
~g  +  C  m2  =  0 

17Tii  771 1  —  1  =  0 

17711  7712  =  0 

17712  7711  =  0 
17712  7712  —  1  =  0 

771 1  17711  +  7712  77712  —  1=0 


where  A,  B  and  C  are  known  and  the  lower  case  letters  a,  6,  c,  e,  /,  g ,  mi,  m2,  imi  and  im2  are  unknown.  We  want  to 
solve  (FAC)  for  the  unknowns  We  use  this  as  an  illustration  because  it  corresponds  to  the  well  known  problem  of  factoring 
a  system  [A,  P,  C,  1]  “minimally”  into  the  product  of  two  systems  (a,  6,  c,  1)  and  (e,  /,  g ,  1).  The  unknown  matrix  (mi  m2) 
corresponds  to  an  isomorphism  between  the  statespace  of  the  (unknown)  product  system  and  the  statespace  of  the  (known) 
original  system.  A  well  known  theorem  says  factoring  is  possible  iff  there  exist  complementary  projections  Pi,  P2  satisfying 

Pi  (A  -  BC)Pi  =  (A  -  BC)Pi  and  P2  AP2  =  AP2  .  (1) 

We  now  apply  a  strategy  to  see  how  one  might  discover  this  theorem.  The  input  is  the  equations  (FAC),  together  with 
declaration  of  A,  P,  C  as  knowns  and  the  remaining  variables  as  unknowns.  Here  is  the  spreadsheet  which  the  computer 
generates: 
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ml**iml**B**C**B**C**ml->A~2**ml  -  A**B**C**ml  - 
SPREADSHEET  B**C**A**ml  +  B**C**B**C**ml  -  ml**iml**A~2**ml  + 

ml**iml**A**B**C**ml  +  ml**iml**B**C**A**ml 


THE  ALGORITHM  HAS  SOLVED  FOR:  {c,g,a,b,e,f > 


The  corresponding  rules  are  the  following: 
a->iml**A**ml  b->iml**B  c->C**ml 

e->im2**A**m2  f->im2**B  g->C**m2 

=======  UNDIGESTED  RELATIONS  APPEAR  BELOW  ===== 

THE  FOLLOWING  VARIABLES  HAVE  NOT  BEEN  SOLVED 
FOR:  {iml , im2 ,ml  ,m2} 

The  expressions  with  unknown  variables 
{iml, ml}  and  knowns  {A,B,C> 


iml**ml->l 

ml**iml**B**C**ml->  -A**ml  +  <=== 

B**c**ml  +  ml**iml**A**ml 


The  expressions  with  unknown  variables  {iml, m2} 
and  knowns  {A} 


iral**m2->0 

iral**A**m2->0 

The  expressions  with  unknown  variables  {im2,ml} 
and  knowns  {A,B,C} 


im2**ml->0 

im2**B**C**ml->im2**A**ml 

The  expressions  with  unknown  variables  {im2,m2} 
and  knowns  {A} 


im2**m2->l 

m2**im2**A**m2->A**m2  <=== 

The  expressions  with  unknown  variables 
{iml , im2 ,ml ,m2}  and  knowns  {B}  * 


m2**im2**B->  B  -  ml**iml**B 
m2**im2->l  -  ml**iml 


The  unknowns  a,6,c,e,/  and  g  are  solved  for.  There  are  no  equations  in  1  unknown.  There  are  4  categories  of  equations 
in  2  unknowns.  A  user  must  observe  that  the  equations  which  I  marked  with  <===  each  transform  to  equations  in  one 
variable  Pi  (respectively  Pi)  when  one  makes  the  assignments: 

Pi  =  mi  imi  and  Pi  —  mi  imi .  (2) 

Run  GBA  again  with  (2)  added  and  Pi,  Pi  declared  known.  The  resulting  spreadsheet  is  much  like  the  one  above  but 
has  the  added  piece 

The  expressions  with  unknown  variables  {  }  and  knowns  {A,B,C,P1 ,P2} 


P2->1  +  -PI  pi**a**P1->P1**A  Pl**B**C**Pl->  -A**P1+P1**A+B**C**P1  P1~2->P1 


P1**A~2**P1->P1**A~2  <===  REDUNDANT 

P1**B**C**B**C**P1->A~2**P1  -  P1**A~2  -  A**B**C**P1  -  B**C**A**P1  + 

B**C**B**0*P1  +  pi**A**B**C**Pl  +  pi**B**C**A**Pl  <“*  REDUNDANT 

You  see  the  first  line  is  the  condition  (1)  of  the  classical  theorem  (plus  stating  that  Pi ,  P2  are  complementary  projections) 
which  immediately  proves  one  side  of  the  theorem.  It  takes  one  more  GBA  run  to  remove  the  redundant  equations  and 
thereby  prove  the  converse  direction. 
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